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Abstract We report on a pilot study on identifying metal-poor stars pre-enriched by 
Pair-Instability Supernovae (PISNe). Very massive, first generation (Population III) stars 
(140 Mq < M < 260 Mq) end their lives as PISNe, which have been predicted by the- 
ories, but no relics of PISNe have been observed yet. Among the distinct characteristics 
of the yields of PISN, as predicted by theoretical calculations, are a strong odd-even ef- 
fect, and a strong overabundance of Ca with respect to iron and the Solar ratio. We use 
the latter characteristic to identify metal-poor stars in the Galactic halo that have been 
pre-enriched by PISN, by comparing metallicites derived from strong, co-added Fe lines 
detected in low-resolution (i.e., R = A/AA ~ 2000) spectra of the Sloan Digital Sky 
Survey (SDSS), with metallicities determined by the SDSS Stellar Parameters Pipeline 
(SSPP). The latter are based on the strength of the Ca II K line and assumptions on the 
Ca/Fe abundance ratio. Stars are selected as candidates if their metallicity derived from 
Fe lines is significantly lower than the SSPP metallicities. In a sample of 12,300 stars for 
which SDSS spectroscopy is available, we have identified 18 candidate stars. Higher res- 
olution and signal-to-noise ratio spectra of these candidates are being obtained with the 
Very Large Telescope of the European Southern Observatory and the XSHOOTER spec- 
trograph, to determine their abundance patterns, and to verify our selection method. We 
plan to apply our method to the data base of several million stellar spectra to be acquired 
with the LAMOST telescope in the next five years. 

Key words: stars: Population II, Population III - stars: supernovae (PISNe) - stars: abun- 
dance - stars: chemically peculiar - method: data analysis - techniques: spectroscopic - 
instrumentation: spectrographs 

1 INTRODUCTION 



The first generation of stars (herafter Population III stars) in the Universe formed several hundred million 
years after the Big Bang. If such stars have survived until today, they can provide us via their abundance 
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patterns with information on the chemical composition of the early Universe. Because the formation of 
Pop. Ill stars occurred in an environment absent of an efficient cooling agent (such as metals and/or dust, 
which are present at later times), they are believed to have been much more massive with charac t eristic 
masses of >=10 M^, up to possibly several hundred solar masses (see, e.g. jBromm et al.ll 19991 



Abel et al.ll2000HKarlsson et alj|201ll). and can re flect the yields of the Big Bang (e.g jAbel et al 



2009; 



2002 



Bromm et al.ll2002l Nakamura & Umemural2002[). This is also suggested by simulations of the collapse 



of primordial molecular clouds ( Ostriker & Gnedinll 19961) . However, some studies have suggested that 
lower mass stars could form within the gas clust er, due to the fragmentation of the tur bulent primodial 
gas in the minihalos (see e.g., IClark et al lEoTTaHbh or due to the protostellar feedback dHosokawa et al.l 
1201 ll) . 

Depending on the mass of the progenit or star, different phy sical mechanisms lead to the explosion 
of Pop. Ill stars at the end of their lives (see lBromm et al.ll2003l) . Their UV radiatio n contributes to the 
re-ionization of the Universe, and they enrich their surroundings with metals (e.g. , iGnedin & Ostrikerl 
1 19971 iFerrara et al.ll2000t iMadau et al.lbOQll iMori et al.ll2002T: lYoshida et al.1l2004 . Very massive first 
stars (140 Mq < M < 260 Mq) will ex plode as Pair-Instability Supernovae (PISNe), leading to com- 
plete disruption of the progenitors (e.g.. iFrver et al.ll2.Q0lt. iHeger & Wooslevl 2002t iLingenfelter et al.l 
120031) . enriching the interstellar medium (ISM,) with heavy elements. Due to their short lives (< 3 My), 
PISNe are likely the first objects in the Universe that have enriched the interstellar medium (ISM). In 
the progenitor stars of PISNe, after central helium burning, the temperature and density in the core sup- 
port pair-instability, and then a partial collapse happens. The collapse proceeds to efficiently compress 
the stars's core, until fast implosive oxygen and silicon burning take place due to the overpressure, re- 
leasing enough energy to r evert the collapse, and the star is complete ly disrupted with no black hole or 
other re mnants left behind (Heger & Woosley|2002l ; lFrver et al1l200lh . The recently observed object SN 
2007bi dGal-Yam et al.ll2009l) is hypothesized to have been a PISN. 

The fraction of metal-poor stars that have formed directly from material that was enriched by the 
yields of PISNe is predicted to be very small dKarlsson et alj|2008h and this together with an obser- 
vational bias may e xplain why no metal-poor star displaying distinct PISNe signatures have been ob- 
served so far (e.g ., IChristlieb et al.l 120021 : ICavrel eta! 1 12004 ICohen et al.ll2004 IBarklem et all 120051: 
iFrebel et"aT1l2005l) . The majority of Pop. II stars with a dominant contribution from PISNe are predicted 
to have Ca abundances in excess of [Ca/H] — — 2.6, which implies th at a significant fraction of the 
PISN-dominant stars may have escaped detection dKarlsson et al.|[20Q8b . In the past, wide-angle spec- 
troscopic surveys for metal-poor stars have used the Ca K line as a metallicity indicator, relying on the 
fact that the v ast majority of the metal-poor stars follow a well-defined trend of [Ca/Fe] as a function 
of [Fe/H] (see lBeers et al.ll 1999b . Due to the high [Ca/Fe] produced in PISNe, the metallicity of Pop. II 
stars pre-enriched by them would therefore be overestimated by about a factor of 40 or more, and hence 
they would not be selected as metal-poor candidates in these surveys. 

The lower side of Fig. 4 o f lKarlsson etaild2008l) shows a simulation of the predicted distribution of 
stars for which > 90 % of the total atmospheric Ca abundance is synthesized in PISNe in the [Fe/Ca]- 
[Ca/H] plane, assuming the number of PISNe occupied 10 % of the primordial stellar population. The 
observation data cover the entire metallicity regime of their model, well reproducin g the small scatter in 
observations (see, e.g. JCavrel et alj|2004llArnone et al.ll20Q5b IBarklem et ai1 l2005). Among the several 
hundred well-studied stars falling in the high number density regions in the abundance ratio diagram, 
none shows the abundance signatures of PISNe. Note that still below [Ca/Fe] ~ —1, there is a negligible 
fraction of PISN-enriched stars. These are the stars that we are targeting with our selection method. 

Non-rotating stars exploding as PISNe differ from core-collapse SNe in that they exhibit a neutron 
excess in their interiors. Consequently, one of the characteristic chemical signatures of PISNe is a strong 
odd-even effect, which means particularly low abundance ratios of odd-Z elements to even-Z elements. 
Furthermore, the production of both rapid and slow neutron-capture elements are absent due to the lack 
of excess neutrons in addition to shorter expansion timescales during the explosion. Other nucleosyn- 
thetic signatures of PISNe are [ Mg,Si/Na,Al] > 0; fS i,S/C1 ~ 1-1.5; [Zn/Fe] < dHeger & Wooslevl 
120021; lUmeda & Nomotoll20oi iKarlsson et alj|201 lb . Hence metal-poor stars pre-enriched by PISNe 
can be identified by these chemical signatures. 
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In the Sloan Digital Sky Survey (SDSS), candidate metal-poor stars are selected by means of ugriz 
photometry; therefore, the above-mentioned bias against PISNe pre-enriched stars is not present in the 
sample of stars for which low-resolution (i.e., R = A/AA ~ 2000) spectroscopy was obtained. 

The seventh data release of SDSS (DR7) contai ns 460,000 stars; for about half of them, low- 
resolution spectra have been acquired. According to iKarlsson et all d2008l) . the number fraction of 
second-generation stars below [Ca/H] = —2.0 with a dominant (i.e.> 90 %) contribution from PISNe 
is 10~ 4 to 5 x 10~ 4 . About 10 % of these PISN-dominated sta rs are predicted to have [Fe/Ca] ~ — 1. 
Therefore, according to the current model dKarlsson et alj|2008l) . at least two stars and up to 9-16 stars 
with dominant pre-enrichment by PISNe could be found in DR7. 

After verifying our method by this pilot study, we will also apply the method to DR8 and fu- 
ture SDSS data releases, as well as to the several million stellar spectra to be obtained in the course 
of the LAMOST (The Large Sky Area Multi-Object Fiber Spectroscopic Telescope) Galactic survey 
(IZhao et al.ll2012h . which will provide more candidates. If we can find any stars pre-enriched by PISNe, 
this would provide observational evidences for their existence; otherwise an upper limit for their space 
density can be derived. 

In this paper we report on the candidate selection applied to SDSS DR7. In section [2] we give a 
brief description of the star sample from SDSS DR7; Section|3]introduces the metallicities from SSPP; 
the data analysis process and candidates selection are described in section |4] and [5] respectively; Finally 
a conclusion is given in section[6] 

2 THE SDSS SAMPLE 

The observations forming the basis of SDSS-DR7 were carried out during an eight-year period with 
a dedicated 2.5-m telescope located at Apache Point Observatory in Southern New Mexico. DR7 in- 
cludes phtometric data of 460,000 stars and over 300,000 stellar spectra covering more than 8200 square 
degrees of the sky. The wavelength coverage is 3800-9200 A, with a resolving power of R = 1800- 
2200 and S/N > 4 per pixel at g = 20.2. We extracted 224,080 stars from the table star joining to 
sppParams in SDSS-DR7 using CasJobs, which is the online workbench for extracting large scien- 
tific catalogs. The number reduced to 17,623 when narrowing the temperature range to 4500-6500 K, 
and requiring the average signal-to-noise ratio (S/N) per pixel of more than 20 in the wavelength range 
4000-8000 A. We adopted this temperature range by inspection of the feasibility of our method on dif- 
ferent spectra. After the data analysis, about 5,000 spectra were rejected due to lack of flux points in 
the wavelength range which contains the iron lines we use. The final sample to which we applied our 
selection criteria hence consists of 12,304 stars. 

When this work was under preparation, the SDSS Data Release 8 (DR8) was published. Besides 
adding more spectra of stars to the data release, the final adopted values of T ff and logg are only 
slightly different from DR7, while the [Fe/H] estimates are in general improved, especially at the low- 
metallicity (< —3.0) and high-metallicity (> 0.0) ends, due to a re-calibration of the NGS1 and NGS 2 
parameter estimation approaches (see TableQ]), and changes in S/N and (g— r)o dSmolinski et al.l201 lb . 
Therefore, the primary sample selection based on T e s and log g is still valid, but we consider the [Fe/H] 
from DR8 instead of DR7. The [Fe/H] from DR7 and DR8 are both listed in Tabled 

3 METALLICITIES DERIVED BY THE SSPP 

SSPP is the stellar parameter pipeline developed for the Sloan Extension for Galactic Exploration 
and Understanding (SEGUE), which is a part of SDSS-II and SDSS-III. This pipeline works on 
low-resolution spectroscopy and ugriz photometry obtained during the course of SDSS-I, SDSS- 
II/SEGUE-1 and SDSS-III/SEGUE-2, making use of multiple techniques to estimate fundamental stel- 
lar atmospheres parameters (including effective temper ature, surface gr av ity, and metallicity), along 
with determining rad ial velocities for AFGK-type stars dLee et al.ll2008allbl; lAllende Prieto et alj|2008l; 
ISmolinski et alj201 ll) . Different approaches were applied for each parameter, as summarized in TableQ] 
For [Fe/H], 12 different methods were adopted in SSPP. 
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Table 1 Metallicity estimation methods used by the SSPP 



VALUE — METHOD DESCRIPTION 



fehl 


NGS2 


X 2 minimization technique using Kurucz NGS a-enhanced grid 1 


feh2 


NGS1 


X 2 minimization technique using Kurucz NGS non-a-enhanced grid 1 


feh3 


ANNSR 


Neural network trained on synthetic spectra 


feh4 


ANNRR 


Neural network trained on real spectra 2 


feh5 


CallKl 


Spectral fitting in the range A = 3850-4250 A, using the NGS 1 grid 1 


feh6 


CaIIK2 


Ca II K line strength combined with a broadband colour 3 


feh7 


CaIIK3 


Neural network approach using the Ca K line index K24 and g — r 1 


feh8 


ACF 


Autocorrelation function proportional to the frequency and strength of weak metallic lines 3 


feh9 


CallT 


Strength of the Ca II triplet lines and B — V; neural network approach 4 


fehlO 


WBG 


Equivalent width of the of Ca II K line and comparison of weaker metallic lines 5 


fehll 


k24 


matching of synthetic spectra using k24 grid of flux with the observed 






flux in the region A = 4400-5500 A, including g — r 6 


fehl 2 


kil3 


matching synthetic spectra using ki 13 grid of flux without colour with 






the observd flux in the region A = 4400-5500 A 1 


feha 


ADOP 


Adopted [Fe/H] value, combination of good estimates 



1 I Lee et al.H2008al 

2 Re Fiorentin et al.|[ 2007 

3 Beers etal.ll 1999 

4 ICena rro et al. 20013 151 

5 I Wilhelm et all 19991 

6 Ullende Prieto et alEwB 

Each of these methods only w orks for a certain colour range and S/N level. More detailed infor- 
mation can be found in Table 5 of ILee et"all d2008al) . Among these approaches, feh 6 and feh7 use the 
strength of Ca II K line in combination with a broadband colour to estimate [Fe/H] dBeers et al l 19991: 
ILee et alJl2008al) . feh9 was obtained by measuring the integrated stren gth of Ca II triple feat ures, which 
are also sensitive to log g, along with the de-reddened B — V colour dCenarro et aLlfeOOhabh . feh5 ap- 
plied the NGS 1 grid to a short wavelength window covering Ca II K, B aimer lines, and a Ca I line, to 
estimate the metallicity, temperature and surface gravity, respectively dLee et alj|2008al) . fehlO uses a 
combination o f the equivalent widt h of Ca II K line and a comparison of weaker metallic lines to syn- 
thetic spectra dWilhelm et al.|[l999h . The adopted feha containing the information of Ca abundance as 
well as other weak metal lines, is much closer to the theoretical definition of metallicity. A flag value 
is assigned to the result from each method if the star does not satisfy the range o f validity for this 
method, or 1 if it is de fined as an outlier after matching with the synthetic spectra (See lLee et alj2008al : 
ISmolinskietai1l201ll) . 

In order to make the comparison more sensitive to the deviation between the abundance of calcium 
and iron, we use the averaged value of feh6 and fehl, if they are both valid, as indicated by a flag value 
of 2; or the valid one of them. Otherwise, one of feh9, fehlO and fehh would be used, depending on 
the given priority, which is based on the correlation degree to Ca abundance. In the case that none of 
them passed the synthetic spectra matching check, the adopted feha was used. 

4 METALLICITY FROM COADDED IRON LINES 

The metallicity Z of a star refers to the proportion of its matter made up of chemical elements heavier 
than helium. Iron is often used as the reference element for metallicity, since it is the end product of 
exothermic nucleosynthesis in stars, and because numerous iron lines are present in the optical wave- 
length region, so that its abundance can easily be measured in high-resolution spectra of cool stars. 
However, if the metallicity is to be determined from low resolution spectra, one has to resort to the 
muc h stronger Ca II H and K lines and assumptions on the relation between [Fe/H] and [Ca/H] (See, 
e.g jBeers et alJll999h . 
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Teff=4500K/logg=4.0 Teff=5500K/logg=2.5 




-30 -20 -10 10 20 30-30 -20 -10 10 20 30 



Relative Wavelength 



Fig. 1 The result of coadding iron lines in synthetic spectra with different [Fe/H] for T e g = 
4500/log g = 4.0 and T e ff = 5500/log g = 2.5, respectively. The line center is shifted to a 
wavevelength of zero in the rest frame. 
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Table 2 The two Fe I lines we coaddcd' 



A [A] 


Species 




log gf 


£[eV] 


4045.82 


Fel 


1.17 


0.280 


1.485 


4383.56 


Fel 


1.01 


0.200 


1.485 



a log(<7/) and lower excitation energy are from Vienna Atomic Line Database (VALD). The equivalent widths 
(EW) are the line strength in the Sun. 
http ://w w w. astro . uni vie . ac . at/ ~ vald/ 



4.1 Selecting iron lines 

We search ed for strong iron l ines detectable in SDSS spectra, using identifications of lines in the Solar 
spectrum dMoore et al.lll966b . Because of the limited quality of the SDSS spectra in terms of resolving 
power and S/N, most of the Fe lines that are very strong in the Sun are difficult to detect in the SDSS 
spectra. After attempts to identify these lines in spectra of stars with different T e g, log g and [Fe/H], the 
lines with equivalent widths below 1 A were rejected. Additionally, because some SDSS spectra have no 
flux data around 3820 A, and also because lines at the extreme blue end with wavelength smaller than 
4000 A are usually blended with other lines or completely hidden in the noise, finally, two relatively 
strong iron lines are left: Fe I 4045.82 A, and Fe I 4383.56 A. Their atomic data are listed in Table [2] 

Fig.[T]displays the coadded iron lines in synthetic spectra of stars with different metallicities, which 
demonstrates that they are strong enough to be detected at [Fe/H]> —4.0 and T e g < 6500K. We hence 
restricted the application of our method to stars in this T e g range. Note that additional information on 
the metallicity is contained in the continuum around the coadded Fe I lines. Our method implicitly 
makes use of this information via the continuum bands that are used to scale the synthetic spectra to the 
observed ones. 

4.2 Normalization 

Continuum normalization of a spectrum is always a delicate task, even more so when the spectra to be 
normalized only have a low resolution, and when the task must be performed automatically for stars 
with a variable range of effective temperatures, which is the case in our application. Because the two 
iron lines we use both lie in the wavelength range of 4000-4500 A, only a short wavelength range in 
the spectrum covering the two iron lines is used. This makes the normalization easier, since a smaller 
wavelength region has to be modeled. We normalize the spectrum in the region A = 3820-5800 A, 
which is a bit larger than the range we use for following anaysis, because we want to include more 
continuum points in the spectrum for a more accurate continuum determination. We adopt the same 
division of the spectrum as that used for two seperate continuum fits in SSPP (Lee et al. 2008), because 
we want to be consistent with the SSPP analysis process as much as possible, which will make the 
comparison of the results more robust. 

In order to determine an approriate continuum fit, some strong lines are removed from the spectrum, 
i.e., Balmer lines, Ca II K lines, Ca I line, Mg lb lines, G-band and iron lines, after which the spectrum 
is iteratively fitted to a ninth-order polynomial. We then rejected the points outside the 3cr range of 
the fitted function after the first fit, which can also remove the cosmic rays. In the following iterations, 
considering the numerous weak lines left in the spectrum, only points within 2a ~ — 0.5cr range of 
the fitted function were kept. The final continuum was obtained after four interations. Fig.|2]shows two 
examples of the normalization for two observed SDSS spectra, among our candidates. 

The synthetic spectra we use were kindly provided by Y. S. Lee (priv. comm.). The wavelength cov- 
erage is 3000-10, 000 A in 0.01 A steps, and they were generated using Kurucz models and the spectrum 
synthesis program TURBOSPECTRUM. The spectra cover a wide range of the stellar parameters T c ff , 
log g and [Fe/H], from 4000 K to 9750 K, 0.0 to 5.0 dex and -4.0 to 0.5 dex, in steps of 250 K, 0.5 dex, 
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Fig. 2 Continuum determination for SDSS spectra of our two candidates: SDSS J0034-0008 
and SDSS Jl 133+4042. 



and 0.5 dex, respectively. The synthetic spectra are processed in exactly the same way as the SDSS 
spectra, in order to reduce systematic errors. 



4.3 Coadding iron lines 

Because the wavelength scale used in SDSS spectra is based on vacuum, the first step before coadding 
the iron lines is to transform the wavelength scale of the original SDSS spectrum to an air-based wave- 
length scale, and after shifting the spectrum to a ze ro-velocity rest frame using the radial velocity es- 
timated from matching ELODI E template spectra jPrugniel & Soubiranll2001[ iMoultaka et alj 120041: 
lAdelman-McCarfhv et ai1l2008l) . the two iron lines are shifted to a wavelength of zero. The radial ve- 
locity from ELODIE is reported to be the best available estimate by previous experiences dLee et al.l 
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l2008ah . Linear interpolation is used to put the spectral regions of the two lines on the same pixel scale, 
before the spectra are coadded. To speed up the analysis process, only ranges of 60 A centered on the 
two iron lines were used. Fig.[3]shows the iron lines detected in the spectra and the strength of the lines 
after coadding, for both observed and synthetic spectra. The synthetic spectra corresponding to the best- 
fit metallicity from our analysis and from SSPP are also plotted for comparison. It is clearly seen that 
the adopted SSPP metallicity is too high for iron lines. It is likely that the SSPP result is dominated by 
other chemical features other than iron, e.g., Ca II, Mg II, or even the molecular carbon bands for some 
cases, which we will discuss more later. 




Fig. 3 Line strength before and after coadding two iron lines of an observed spectrum for 
one of our candidates (solid, SDSS J21 18 — 0613) and the synthetic spectra (dashed, T c ff = 
4750K/logg = 4.5/[Fe/H] — — 2.0/0. Odex). The red and green dashed line correspond to the 
best fit synthetic spectrum from our analysis ([Fe/H]=— 2.04) and from SSPP ([Fe/H]=— 0.10), 
respectively. 



4.4 Metallicity determination 

Based on the T ff and logy from SSPP, the corresponding synthetic spectrum was chosen, and then 
it was degraded and interpolated to the same resolving power (R ~ 2000) and wavelength scale as 
SDSS spectrum. Since the resolution does not change much in the wavelength range we use in this 
work, we adopted a constant wavelength step width. Then we used a reduced x 2 criterion to get the 
best-fit metallicity by minimizing the difference between the observed and synthetic coadded flux. We 
used a range of 10 A of the coadded spectrum, which only covers the coadding feature, centered on 
the line, in order to reduce the impact from other lines and noise in the rest part of the spectra. A 
third-order polynomial is adopted to fit all the \ 2 values from different metallicities, and the metallicity 
coresponding to the minimum \ 2 is taken as the best-fit value. 

4.5 Error estimation 

In order to estimate the uncertainty of our method, we analyzed 1069 stars belonging to different globu- 
lar clusters (GCs) and open clusters (OCs). The spectroscopic data were obtained during SEG UE obser- 
vations and they were used to check the validation of the SEGUE target selection and SSPP (Lee et al. 
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l2008bt ISmolinski et alJl201ll) . Only GC members with 4500 K< T cff <6500K and S/N > 20 were 
selected in consistency with our sample, as described in Section [2] Finally, 724 stars were considered 
for eight GCs (M3, NGC 5053, M53, M71, M92, M15, M2, and M13) and five OCs (NGC 6791 M35 
M67, N GC 2420, a nd NG C 2 158). The high-resolutio n metallicity results are from the work of lHarrisI 
d!996l) : iDias etal] d2002l) and iBoesgaard et al.1 d2009l) . By comparing the metallicities from our work 
with the averaged [Fe/H] from high-resolution analysis for each cluster, we estimated the uncertainty of 
our method to be ~ 0.40 dex. 

5 CANDIDATES SELECTION 

Due to the differences in the methods used, together with the weakness of the Fe I lines, the limited qual- 
ity of the SDSS spectra, and uncertainties incured by the automated analysis, our [Fe/H]_x 2 estimated 
directly from Fe I lines expectedly deviates from the SSPP metallicities [Fe/H]_sspp. The average devi- 
ation is 0.29 ± 0.25 dex. As shown in Fig.|4j at the high and low metallicity ends, larger scatters appear 
in the comparison. A linear fit [Fe/H]ggpp = 0.76 x [Fe/H] ^1^0^ — 0.23 is performed, with a stan- 
dard deviation of 0.35 dex, which demostrates that our result is systematically higher at [Fe/H]> —1.0. 
The 3cr outliers based on the linear fit are interesting, despite possible origination from errors. In this 
work, we care about those objects falling outside 3er with likely Ca-enhancement, which might be po- 
tential candidates for PISNe enriched metal-poor stars. By means of visual inspection of the original 
and coadded spectra, we rejected false positives from the candidate sample, i.e., spectra suffering from 
overestimated signal-to-noise ratios or poor continuum determination. Finally, 18 stars were selected as 
Ca-enhanced candidates. They are listed in TableO and shown as red stars in Fig. |4] 

The spectra of outliers with lower SSPP metallicity were also inspected, and we found that the large 
departure is mainly due to low S/N and hence poor continuum determination. 

In addition, the ou tliers are mostly c ooler (T c ft < 5000 K) and metal-poor stars ([Fe/H]< —1.00), 
therefore, according to iLee et all d2008al) the estimated metallcity from Ca lines in this case might be 
underestimated by ~ 0.5 dex, indicating the true departure from our derived results would be even 
larger. 

Although We note that a large fraction of the used [Fe/H] values from the SSPP are the final adopted 
value feha, which contains the information of Ca and other metallic lines, e.g., Mg I, Fe I, Na I, the true 
difference (excluding uncertainty) between the two metallcities can still reflect the deviation of abun- 
dance between Fe lines and other metal lines, especially Ca lines, since some weak Ca lines might be 
used and the final valid metallicity estimators were selected after matching with th e synthetic spectrum 
in the wavelength range of 3850-4250 A and 4500-5500 A dSmolinski et al.ll201 lb . 

We inspected Ca I 4226 A , Ca II H and K lines in the spectra of our candidates, and no apparent 
enhancement of Ca abundance was found on Ca I 4226 A line, while for Ca II H and K lines, the situa- 
tion is complicated due to the saturation and blendings. Furtherly, we are not clear about the sensitivity 
of Ca abundance on the calcium line strength, and whether these stars are influenced by chromospheric 
activities, thus better quality data are needed to measure the accurate abundances and look into these 
questions. In Fig. [5] our metallicity results also indicate that there are larger differences bewteen the 
[Fe/H] values of DR7 and DR8 at the low metallicity tail compared to the differences at higher metal- 
licity. The number of stars with a large offset in the comparison between DR8 and our derived [Fe/H] 
increases when a smaller number of estimators is used, and fewer estimators (< 3) are used at lower 
metallicity ([Fe/H] < -3.0). 

In general, our method to derive [Fe/H] is reliable, as seen from their good correlation with the 
metallicity of different approaches used in SDSS-DR8. We also note that only a few selected outliers 
have valid result from methods using Ca II lines, which weaken the correlation with Ca-enhancement. 
We do not exclude the possibility that the large departure of our [Fe/H] from those derived with the 
SSPP is caused by errors in either or both of the two estimates. Therefore, it is necessary to acquire 
spectra of higher resolution and S/N to verify our selection method, and for possibly identifying via 
their abundance patterns stars pre-enriched by PISNe. To this end, medium-resolution (R ~ 6000- 
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Fig. 4 Comparison between metallicities from SSPP and our method. Those 3<j outliers with S/N > 25 
(red stars) on the upper side show a possible Ca enhancement. Both the upper and lower 3<r outliers are 
highlighted in red dots. 



11, 000) of our candidates are currently being obtained with the XSHOOTER spectrograph attached to 
the Unit Telescope 2 of the Very Large Telescope of the European Southern Observatory. 

6 CONCLUSIONS 

Very massive (140 Mq < M < 260 Mq), first generation stars end their lifes as PISNe, which have 
been predicted by theories, but no strong observational evidence for their existence in the early Universe 
has yet been obtained. The theoretically predicted chemical fingerprints of PISNe are high overabun- 
dances of Ca with respect to Fe and the ratio of these elements in the Sun, and a strong odd-even effect. 
We carried out a pilot study on identifying metal-poor stars with strong Ca-enhancement, by comparing 
the metallicities determined with the SDSS SSPP containing information on the Ca abundance (i.e., 
those estimators which involve the Ca II H and K lines or the Ca II infrared triplet) with metallicities 
we derived directly from iron lines. We used a coadding technique to increase the S/N of the SDSS 
spectra. A \ 2 minimization method was applied to search the best matched metallicity from a series of 
synthetic spectra with known atmosphere parameters. Our derived metallicities correlate well with those 
of the SSPP, with a linear fit [Fe/H]sspp = 0.76x [Fe/H] t hiswork — 0.23. By applying a 3er criterion and 
visual inspection, 18 candidates were selected. Once the nature of these candidates is determined and 
our selection method is verified, we will apply our method to the low-resolution stellar spectra obtained 
with the LAMOST telescope. 
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Fig. 5 Left: Metallicity deviation of DR7 from DR8, distributing with our derived metallicity. 
Right: Metallicity deviation of DR8 from our method, distributing with the number used to 
give the final metallicity in SSPP DR8. Black dots: the whole sample; red circles: 3d outliers 
with higher Ca; red filled circles: selected candidates. 



Table 3 Candidates of Metal-poor Stars Pre-enriched by PISNe. 



Star 


Plate 


MJD 


Fiberid 


S/N 


Tcff 


log 9 


[Fe/H] 


[Fe/H] a 


[Fe/H] 


Note 6 












(K) 




DR7 


DR8(Ca) 


this work 




SDSS J081157. 13+143533.0 


2270 


53714 


253 


52 


4750 


4.0 


-0.95 


-1 


20 


-3.52 


adop, 1 


SDSS J015315.91-084216.8 


0665 


52168 


489 


35 


4822 


3.6 


-1.00 


-0 


96 


-2.98 


adop, 1 


SDSS J21 1851.04-061319.6 


0639 


52146 


538 


47 


4795 


4.4 


-0.13 


-0 


10 


-2.04 


adop, 2 


SDSS J030209.90-062704.0 


0458 


51929 


365 


33 


4657 


4.2 


-0.62 


-0 


65 


-2.43 


adop, 1 


SDSS J121 144.01+01 1951.3 


0517 


52024 


003 


65 


4798 


4.5 


-0.19 


-0 


17 


-1.79 


adop, 2 


SDSS J205923.08-065208.2 


0636 


52176 


064 


42 


5298 


4.0 


-0.91 


-0 


60 


-2.48 


wbg, 3 


SDSS J080805.07+183352.0 


1923 


53319 


595 


31 


4830 


4.8 


-1.19 


-0 


99 


-2.76 


caiik2, 8 


SDSS J 1 1 223 8 . 1 8+4345 55.3 


1365 


53062 


144 


26 


4808 


4.5 


-1.11 


-0 


77 


-2.65 


caiit, 6 


SDSS J164435.94+203659.1 


1569 


53168 


232 


36 


4665 


4.6 


-0.32 


-0 


32 


-1.66 


caiit, 6 


SDSS J075800.22+435257.9 


0437 


51869 


374 


60 


4588 


4.9 


-0.13 


-0 


02 


-1.37 


adop, 1 


SDSS Jl 11218.53+262543.8 


2212 


53789 


636 


43 


4624 


5.0 


-0.09 





00 


-1.31 


adop, 3 


SDSS Jl 13344.77+404238.4 


1443 


53055 


090 


36 


4549 


4.7 


-0.61 


-0 


19 


-1.83 


caiikl, 5 


SDSS J155219.79-001022.7 


0342 


51691 


038 


35 


4662 


3.9 


-1.70 


-0 


99 


-2.78 


caiik3, 2 


SDSS J104459.32+360554.7 


2090 


53463 


273 


33 


5877 


4.3 


-0.95 


-0 


69 


-2.66 


adop, 3 


SDSS J095423. 10+240410.8 


2341 


53738 


019 


27 


5004 


3.6 


-1.93 


-1 


16 


-3.12 


adop, 1 


SDSS J230130.60+001920.9 


0380 


51792 


499 


30 


4798 


4.5 


-1.08 


-1 


17 


-2.96 


adop, 2 


SDSS J153642.56+210747.4 


2166 


54232 


266 


31 


5489 


4.2 


-1.02 


-1 


30 


-2.85 


adop, 1 


SDSS J0034 16. 18-000849.9 


0392 


51793 


224 


27 


5462 


4.1 


-0.75 


-0 


61 


-1.91 


wbg, 3 



a : [Fe/H] of DR8(Ca) refers to the Ca-related metallicity from SSPP(DR8) used in this work. 

b : This column list the method used to determine [Fe/H]-DR8(Ca) and the number of valid methods in SSPP. 
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